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The contributions of confining as well as nonconfining nonperturbative self- 
interactions of stochastic background fields to the shear and bulk viscosities of the 



o 

gluon plasma in SU(3) Yang-Mills theory are calculated. The nonconfining self- 



interactions change (specifically, diminish) the values of the shear and bulk viscosi- 
ties by 15%, that is close to the 17% which the strength of the nonconfining self- 
interactions amounts of the full strength of nonperturbative self-interactions. The 

i ■ 

Q_i' ratios to the entropy density of the obtained nonperturbative contributions to the 

<D : 

shear and bulk viscosities are compared with the results of perturbation theory and 
the predictions of TV = 4 SYM. 
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"tf- ■ I. INTRODUCTION 

<N 

Ultrarelativistic nucleus- nucleus collisions performed at RHIC indicate that, in the vicin- 
ity of the deconfinement critical temperature, the quark-gluon plasma can behave more like 
liquid than a dilute gas of quarks and gluons. One indication of this kind is the experi- 
mentally found common velocity of different species of particles, which are emitted from tie 
expanding fireball of the quark-gluon plasma. The corresponding experimental results 
can be successfully described by the relativistic hydrodynamics of an almost non-viscous 
liquid [3j. The energy-momentum tensor of an ideal (i.e. absolutely non- viscous) liquid has 
the form 4j 0^ = —pg pv + Tsu p u v , where s = s(T) is the entropy density, and u p is the 
velocity of energy transport. The principal deviation from the ideality is described by a 
correcting term AO^, which depends on the derivatives of the velocity linearly: 



A6^ = i] T ■ (A p u u + A v u p ) + - Cr^j H pu d p u p , 

where H pu = u p u v — g pu , A p = d p — u p u v d v . The temperature-dependent coefficients, t]t 
and Ct, are called respectively the shear and the bulk viscosities. Together, they are called 
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first-order transport coefficients of the liquid. 

In hydrodynamics 4j, the shear viscosity characterizes a change in shape of a fixed vol- 
ume of the liquid, whereas the bulk viscosity characterizes a change in volume of the liquid 
of a fixed shape. In particular, the bulk viscosity enters the Navier-Stokes' equation through 
the term (£r + ^Vt) graddivv, being therefore relevant only when div v 7^ 0, i.e. for com- 
pressible liquids. The shear viscosity enters the Navier-Stokes' equation foremost through 
the term 77^ Av, representing the ability of particles to transport momentum. Experimental 
data for water, helium, and nitrogen, analyzed recently in Ref . [5] , show the minimum of the 
shear- viscosity to the entropy-density ratio (j]t/s) and the maximum of the bulk- viscosity 
to the entropy-density ratio (Ct/s) in the vicinity of the liquid-gas phase-transition critical 
temperature. Given different types of liquid-gas phase transitions and different types of 
molecules for these three substances, one can guess that the temperature-behaviors of t]t/s 
and Ct/s are universal, and thus can be qualitatively the same also for the quark-gluon 
plasma. 

Furthermore, it is worth noticing that the (r^/^-ratio has been calculated in M = 4 
SYM theory |6|], where it equals to the temperature-independent constant 1/(4tt). [Note 
that, since M = 4 SYM theory is a conformal field theory, its deconfmement phase can start 
only right at T = 0, for which reason one should not expect any temperature-dependence 
of the (r/r/s) -ratio.] This constant is at least by one order of magnitude smaller than the 
minima of the above-mentioned empirical results for water, helium, and nitrogen, as well as 
the values of the (?7<r/s)-ratio in perturbative QCD. Such a smallness of this ratio in M = 4 
SYM has led to the widely known conjecture that the finite-temperature version of M = 4 
SYM provides (a theoretic example of) the most perfect quantum liquid. For example, at the 
temperature of 200 MeV, the length of a mean free path of a parton traversing such a liquid 
is as small as <0.1 fm. However, since M = 4 SYM is a conformal field theory, the bulk 
viscosity (t in this theory is strictly zero. That makes M = 4 SYM different from the real 
QCD, where the effects of non-conformality are essential (e.g. in the so-called interaction 
measure e — 3p) up to the temperatures as high as 2 3 times the deconfmement critical 
temperature T c . 

In a quantum field theory, the contributions of various fields to the viscosities can be 
calculated by means of the Kubo formulae. A Kubo formula is an integral equation, 
which expresses the spectral density pr{oo) of a given viscosity via the corresponding finite- 



temperature two-point Euclidean correlation function of the energy-momentum tensor 
specifically (612(^)012(2/) ) T f° r t ne shear viscosity and (6 AtM (x)0 1/i/ (?/)) T for the bulk vis- 
cosity. In particular, the fact that it is the trace anomaly 8^(1), which is correlated in case 
of the bulk viscosity, explains in terms of the Kubo formula why the bulk viscosity vanishes 
in any conformal field theory Below, we study a more realistic case of SU(3) Yang-Mills 
theory. There, the energy-momentum tensor can receive contributions from the stochastic 
background fields and the so-called valence gluons. The low-energy nonperturbative stochas- 
tic background fields are characterized by the temperature-dependent chromomagnetic gluon 
condensate (g 2 (F^) 2 } T and the correlation length of the chromomagnetic vacuum, 1//It, that 
is the distance at which the two-point correlation function of the chromomagnetic fields ex- 
ponentially falls off jsl, I2I. The valence gluons are higher in energy than the background 



fields, and are confined by these fields at large spatial separations 



lOj . This property of the 



valence gluons makes them different from the ordinary gluons of perturbation theory. Such a 



two-component model of the gluon plasma was recently proven efficient (cf. Refs. [H|, 12]) 
for the description of the pressure and the interaction measure, which had been simulated 
on the lattice in Ref. [13], as well as for the calculation of the radiative energy loss of a 
highly energetic parton traversing the plasma jl4|. 

It therefore looks natural to apply the same two-component model of the gluon plasma to 
a derivation of the shear and bulk viscosities by means of the Kubo formulae. The first step 
in this direction was made in Ref. jl5|, where a contribution of confining self-interactions 
of the stochastic ba ckg round fields to the shear viscosity was found. Besides those, lattice 



simulations [9|, LL6|, LL7J point to the existence of also nonconfining nonperturbative self- 
interactions of the background fields, whose strength amounts to 17% of the full strength 
of nonperturbative self-interactions. In this paper, we calculate the contribution to the 
shear viscosity produced by nonconfining nonperturbative self-interactions of the background 
fields, and show that it diminishes the contribution of confining self-interactions by 15% (that 
is rather close to the above-mentioned 17%). We furthermore calculate for the first time the 
contributions of both confining and nonconfining nonperturbative self-interactions to the 
bulk viscosity of the gluon plasma. 

Note that the Kubo formulae lead to the following temperature dependence of the con- 
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tributions of the background fields to the viscosities 



id: 



2 



rj T oc Ct oc f L . 

At temperatures above that of the dimensional reduction, T > T* (where T* lies somewhere 
between T c and 2T C |l3|), one has i]t oc (t oc {Ot^) 3 , whereas the temperature dependence 
of the entropy density at T > 2T C is s oc T 3 . Therefore, at such rather high temperatures, 
one expects that the ratios to the entropy density of the contributions of the background 
fields to the two viscosities exhibit the same temperature behavior: 

^cx^oc^ at T>2T C . 
s s 

In this paper, we prove that this is indeed the case, and explicitly calculate numerical 
coefficients in the formulae above. The contributions of valence gluons to the shear and bulk 
viscosities will be addressed in the future publications. 

The paper is organized as follows. In the next section, we first review our approach, and 
then use it to calculate the contribution of confining self-interactions of the background fields 
to the bulk viscosity of the gluon plasma, Ct- in section III, we find contributions to both 
t]t and (t produced by nonconfining nonperturbative self-interactions of the background 
fields. In section IV, we perform a numerical evaluation of the obtained contributions to the 
viscosities, and compare their ratios to the entropy density with those known from other 
approaches. In section V, we summarize the results of the paper. 



II. CONTRIBUTIONS OF CONFINING SELF-INTERACTIONS OF THE 

BACKGROUND FIELDS 



We start with recollecting the approach to the calculation of the shear viscosity 7/t, 
suggested in Ref. jl5|, and generalizing it to the calculation of the bulk viscosity (t- The 
two viscosities can be defined through their spectral densities, p^\co) and pP(oo), as Q 



Vt 



7T- 



dp, 



duj 



(6) 



and 



7r dp T 



uj=Q 



9 du 



(1) 



uj=0 



Here, the superscripts (s) and (b) stand respectively for "shear" and "bulk". Both spectral 
densities can be obtained from the integral equation, called Kubo formula: 

cosh [tv (2:4 — f )] 



dupP' {b \co) 



sinh(w/3/2) 



d 3 x U^ sUb) (x,x 4 -f3n). 



(2) 
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In this equation, Uj! and are the following correlation functions of the Yang-Mills 



energy-momentum tensor 



4 S W 4 - M = <0i2(O)e 12 (x,x 4 - /3n)) T , where 6 12 = g 2 F^F 2fl , 



4 b) (x,x 4 - 0n) = <e w (O)0 w (x,x 4 - 0n)) T , where 9 W = ^-(F« u ) 2 . (3) 

^9 



In Ref. 



15), Eq. 02]) for p^\u) was solved by means of the Fourier transform. Owing to the 
uniformity of Eq. ([2]), this approach equally applies to p^\u>). Its main idea is to use for the 
nonperturbative parts of the zero-temperature correlation functions (x) the following 

ansatz, exponentially falling off with distance: 

umx)=NL .m (G ^.!^m. (4) 

Here (G 2 ) = (g 2 (F^ u ) 2 ^ is the gluon condensate, K 2 - a is the MacDonald function, and 
a G (0, oo) is some parameter. Assuming parametrization fll]), one gets at T > T c the 
following Fourier transform of Eq. ([2]) : 

poo T,/r2a—i 

where T(a) is the Gamma-function. Note that, while the gluon condensate and the corre- 
lation length 1/M become temperature-dependent, the overall coefficients Na^'^ are deter- 
mined entirely by the corresponding zero-temperature correlation functions (x) . Next, 
one solves Eq. fl5]) by imposing for p^ (uj) a Lorentzian-type ansatz 

p^ {b \to) = d° m r . (6) 

For any a G (0, oo), it provides convergence of the (^-integration in Eq. ([5]), and guarantees 
that both sides of Eq. ([5]) have the same large- \k\ behavior. Furthermore, it turns out that 
only for a single value of a, namely a = 1/2, the spectral densities, once sought in the form 
of Eq. ([6]), appear /c-independent 15]. [As can be seen from Eq. Q, at a = 1/2, the spectral 
densities take the purely Lorentzian form considered in Refs. 7|, LL8|.] The corresponding 
temperature-dependent functions C^'^ can then be found, and read 
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Thus, to find the viscosities, it remains to determine the coefficients N^ a l . To this 
end, one evaluates the correlation functions U^ Q (x) = (©i2(0)©i2(^)) and U^L (x) = 
{® titi(Q)®vu(%)) y i & the so-called Gaussian-dominance hypothesis. This hypothesis, sup- 
ported by the lattice simulations {9], Q, states that the connected four-point cor- 
relation function of gluonic field strengths can be neglected compared to pairwise prod- 
ucts of the two-point correlation functions. For the correlation function U^ (x) = 
(^ 4 F 1 a M (0)F 2 a M (0)F 1 6 I/ (x)F 2 ^(x)), this approximation yields 

U%U(x) * <^(0)F£(0)) {9 2 Fl{x)Fl{x)) + (^(0)^)) (g 2 F 2 ^(0)Fl(x)) + 

+ (g'F^FKx)) (g 2 F^(0)Fl(x)) . (8) 

The sought nonperturbative contributions to the viscosities, produced by the background 
fields, can now be parametrized by means of the stochastic vacuum model 8]. In this 
section, we consider only the contribution produced by confining self-interactions of the 
background fields, while the contribution of nonconfining nonperturbative self-interactions 
will be considered in the next section. Confining self-interactions of the background fields 
enter the two-point correlation function of gluonic field strengths as Jg[ [ijj] 

(C 2 ) t) ab 
(g 2 F; u (x)F b Xp (0)) = ^ ■ (5, x 6 up - <WW ■ ■ D(x). (9) 

For the rest of this section, we set N c = 3. At large \x\, the dimensionless function D(x) 
falls off as e~ M ' x ', where \ is the so-called vacuum correlation length js, 9, 16| (for a review 



sec 



17]). The compatibility of Eqs. (jSJ) and (j9J) with Eq. (j3J), at a = 1/2, is achieved by 



having M = 2\x and 



Plugging Eqs. and f lTUj) into Eq. (jSJ), one obtains the desired coefficient N^ s l in terms of 
the constant A: 

N\% = — . (11) 

1/2 57 g V ) 

The constant A can be fixed by means of the expression for the string tension in the funda- 
mental representation 

0t = H4 / d2xD{x) > (12) 
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which yields 



A 



1.05. 



(13) 



J" °° dz ■ zV* ■ jK^Az) 

Finally, substituting Eq. (fill into Eq. (JJj), and using Eqs. (CQ) and ([6]), one accomplishes the 



calculation of the shear viscosity 



15|: 



7T 



1]T 



V 2 A 2 (G 2 ) 



(14) 



4608^ /4 

A contribution of confining self-interactions of the background fields to the bulk viscosity 
can be found in a similar way. It amounts to calculating the coefficient Ny 2 in Eq. (jlj). 
That can be done by using in Eq. (j3J) the one-loop Yang-Mills ^-function, in which case 

2 



2g 



■-^g 2 , and U?Jx) ~ 



11 

32^ 



x 



(15) 



Furthermore, by means of the Gaussian-dominance hypothesis, the correlation function 



(x) can be written as 



U? =0 (x) 



11 

32^ 



G 2 ) 2 + 2(g 2 F« u (0)K( x )Y 



(16) 



Applying the parametrization of the stochastic vacuum model, Eq. ([9j), we have 



U? =0 (x) 



11 

32^ 



G 



2\2 



1 + — L> 2 (x) 
24 v ; 



(17) 



The renormalized spectral density corresponds to UfL Q {x) with the "1" in the brackets of 
Eq. (117j) subtracted. The function D(x) in the form of Eq. ( ITUl) yields then the coefficient 
which enters Eq. (J4]): 

" 2 24 V32i 2 , 

Plugging this coefficient into Eq. 0, and using Eqs. ([I]) and ([6]), we arrive at the following 
formula for the bulk viscosity: 



A 2 



in 2 (g 2 ) 



- 



17287273 \32J /4 
Its numerical evaluation will be performed in section IV. 



(18) 
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III. CONTRIBUTIONS OF NONCONFINING NONPERTURBATIVE 
SELF-INTERACTIONS OF THE BACKGROUND FIELDS 

The two-point correlation function {g 2 F^ u (0)F^ (x)") parametrizes not only confining self- 
interactions of stochastic background fields, through the function D(x), but also nonconfin- 



ing nonperturbative self-interactions, through the so-called function Di(x) 
two-point correlation function reads [cf. Eq. ([9])] 

5 ab 



19|. The full 



(^(0)^>)> = ^ 



N? 



— ■ {k(^\5 up - 5 lip 5„x)D(x)+ 



+ 



[d^xxdup - x p 5 u x) + d u (x p 5px - x x Sp P )] D^x)}, 



(19) 



where k G [0, 1] is some parameter. 

To see that the functions D(x) and Di(x) indeed parametrize respectively the confin- 
ing and the nonconfining nonperturbative self-interactions, one can substitute parametriza- 
tion (JT9~j) into the Wilson loop [21] (W(C)} = (tr P exp (ig § c dXpT a A a ^j ) expressed through 
the correlation function ( g 2 F"(x) F\ (x 1 )) by means of the non-Abelian Stokes' theorem 
and the cumulant expansion p, ll9|. That yields (cf. the second review in Ref. 8j) 

C 2 (G 2 ) 



(W(C)) = exp 



96(iV c 2 - 1) 
+ (1 — k) <t> dx 



2k 



d(Jpv(x) 



da^(x')D(\x - x'\) + 



dx' 



^i(v^) 



(20) 



<C JC ' J(x-x') 2 

where C 2 is the quadratic Casimir operator of a given representation (i.e. T a T a = C2I). 
From this expression, it is explicitly seen that the function D(x) mediates confining self- 
interactions of the background fields, as described by the double surface integral, while the 
function D\(x) mediates nonconfining nonperturbative self-interactions, as described by the 
double contour integral. These two functions can be viewed as phenomenological propagators 
of the background gluons, which describe the two types of self-interactions. In what follows, 
we set 

D 1 (x) = D(x), 

as suggested by the lattice data js), Q, 17 1. This assumption means that, with the separation 
between two points in Euclidean space, the nonconfining nonperturbative self-interactions 
fall off at the same vacuum correlation length l//i as the confining ones. The two types 
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of self-interactions differ, though, in magnitude, that is taken care of by the parameter k. 
In particular, at k — 1, Eq. (Tl9|) reproduces Eq. ([9]), implying the pure area-law of the 
Wilson loop and the full suppression of the nonconfining nonperturbative self-interactions. 
Instead, the opposite limit, k = 0, describes a nonconfining vacuum, in which the Wilson 
loop respects the pure perimeter-law. The lattice-simulated value of k in SU(3) Yang-Mills 



theory is 



17| k = 0.83 ± 0.03. Below, by setting k, to this realistic value, 



K ~ 0.83, (21) 

we relax the approximation k = 1 adopted in the previous section and in the earlier works 
15] . This way, we account for the contributions to the shear and bulk viscosities produced 
by nonconfining nonperturbative self-interactions of the background fields, whose strength 
thus amounts to 17% of the full strength of nonperturbative self-interactions. 

We start with the bulk viscosity, and recalculate the correlation function U^L^x), 
Eq. ( TT6l) . by using for (g 2 F" u (0) F% (x)} parametrization ( fl9l . We denote the correlation 
function with k / 1 as ^r=o( x )- ^ straightforward calculation yields for it the following 
expression [cf. Eq. ((T 



Mb) 

>T=0\ 



11 

32^ 



G 2 > 2 <|1 



1 

24 



—DD' + 2l 



f\2 



(22) 



where D' = We seek the function D(x) in the form generalizing Eq. (|10J) : 



= A ■ / K (yu|x|), where A K =i = A and / K =i(z) 



'#3/2(2*) 
(2*) 3 / 2 



(23) 



Given that U^ (x), Eq. (j22|) . contains terms to the order (1 — k) 2 , the function f K = f K (fi\x\) 
can be sought as a power series 

/„ = f K=1 + (1 - + (1 - K ) 2 / {2) + O ((1 - «:) 3 ) • 

Furthermore, the Matsubara-mode independence of the spectral density is achieved 
when the correlation function (|2"2"|) (with "1" in the curly brackets subtracted) has the form 
of Eq. (J1J) with a = 1/2. The equation representing this condition, 



fl + 



1 - « 



fkf k + 



if, 



r\2 



Ks/2(Mx\} 
(2/i|x|) 3 / 2 ■ 
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yields the following functions f^(/j,\x\) and f^ 2 \jjL\x\): 



f(i) = _if fP) = — 

J K=ll J 

where again /' = |x|J^t. Explicitly, we find 



hi _ (f K =i) 21 
Jk=1 2f K=1 



7T 1 / 4 e~ z 
= 2567257? >(! + W -mi + ^f+ 

+ (1 - «) ■ 16(1 + 2z)(3 + 6z + Az 2 ) + (1 - k) 2 ■ [9 + Az ■ (9 + z ■ (7 + 4z(l + 2)))] }, (24) 

that reproduces correctly the function 

vr 1 / 4 e~ z 
/^) = ^-^-(l + ^ 1/2 - 

given by Eq. (1231) . Next, much as the coefficient .A, the coefficient A K can be found by means 
of relation (fT2|) . and reads 

J dz-z-f K {z) 

[In particular, using the explicit form of f K=1 (z), Eq. (|23p . one can see that Eq. (|25)) at 
k = 1 reproduces Eq. (II 3 p . ] Plugging into Eq. (j25p the obtained function f K , Eq. §M§ , at 
the lattice-simulated value of k, Eq. (l2Tj) . we find 

A=o.83 ~ 0.97. (26) 
The corresponding bulk viscosity at this value of k, 

,.=0.83 _ (A=q. 83 ) 2 fnV (g 2 ) 2 t 

appears by a factor of 



>^k=0.83 X 



( , - 0.85 (28) 

smaller than the approximate bulk viscosity ££ =1 given by Eq. (1T81) . 

We proceed now to the shear viscosity, and recalculate the correlation function U^l (x) 7 
Eq. (jSJ), using for \g 2 F® (O)i^ parametrization ( |T9|) . Each of the condensates 
<^ 2 F 1 a /J (0)F 2 a ^(0)) and <^ 2 F 1 6 i/ (x)F 2 ^(x)) on the RHS of Eq. QH} can only be propor- 
tional to 5\2, and therefore vanish. Upon the use of parametrization (fT9l) . the terms 
<2 2 F^(0)Ft(:r)) (p a i^(0)Jf* (ar)) and (^(0)^)) (g 2 F^(0)Fl(x)) read 



(^(O)^)) <^ 2 F 2 ;(0)F, 



2u( X 
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(G 2 ) 



2\ 2 



72(iV ( 



D 



-DD' 



:i-k) 



X\X2 



,,y 

\x\ J 



X 



and 



(9 2 F?,(0)Fl(x)) (g 2 F^(0)Fl(x)) 
The sum of Eqs. ((29]) and (j3QJ) yields 



(G 2 > 



2\ 2 



72(iV c 2 - 1) 



X\X 2 

\x\ 



dD 

d\x\ 



dD 

d\x\ 



(29) 



(30) 



4t (a? - 



D 2 + ^—-^-DD' + ^ 



72(iV c 2 



■ x 



i Z Z 

+ X — X-i — X 




(31) 



Furthermore, since the self-interactions of stochastic background fields we are studying are 
nonperturbative, it is legitimate to use in Eq. (13T|) the leading large- \x\ approximation. In this 
approximation, the sum (x\ + x 2 ) can in general be disregarded compared to x 2 . We should 
then compare x 2 with 4(xiX2/\x\) 2 , that is the same as to compare x 2 with 2|xiX2|. This 
comparison obviously leads to the inequality (|a;i| — |a;2|) 2 + x\ + x\ > 0, which means that, 
in the leading large-|x| approximation, the term 4(xiX2/\x\) 2 can be disregarded compared 
to x 2 . Equation ( 13~TT) then takes the form 

(G 2 ) 2 



D 2 + 



-DD' 



'1 



t\2 



72(N 2 — 1) [ 2 ' 8 

Comparing this expression with Eq. (1221) . we conclude that the function D(x) has the same 
form of Eq. ( 1231) as for the bulk viscosity, with the function f K (x) and the coefficient A K 
given by Eqs. (1241 and (125]) . respectively. Thus, we obtain the following expression for the 
shear viscosity with the nonconfining nonperturbative self-interactions of the background 
fields taken into account (at iV c = 3): 



'It 



vr 5 / 2 (A=o. 83 ) 2 (G 2 ) 



(32) 



4608\/2 

It is by 15% [cf. Eq. (128]) ] smaller than Eq. (fl4l) . where only confining self-interactions are 
taken into account [15]. In the next section, we proceed with the numerical evaluation of 
the ratios to the entropy density of the bulk and shear viscosities, respectively Eqs. (]T8]) . 
$27]), and ([32]). 
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IV. NUMERICAL EVALUATION 



For the calculation of the bulk viscosity, we use the one-loop running coupling [cf. 
Eq. flUD] 

-2 

9t 



, , T , , 11N C 
2ooln— , where o 



487T 2 ' 

and, for the (N c = 3)-case under study, A = 0.104T C , T c = 270 MeV [J3j. This coupling 
should be plugged into the function 



f(T) 



1 at T c <T < T», 

9t„ ' 1 * 



(33) 



which defines the temperature-dependent inverse vacuum correlation length and the chro- 
momagnetic gluon condensate as Jj], 14. fis| 



M'/CO, <^ 2 > T = (G 2 ) ■ f\T). 



The corresponding zero-temperature values of these quantities are /i = 894 MeV [9] and 
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(G 2 ) = ■^■Vf^ 2 — 3.55 GeV 4 [19] . where cr f = (440 MeV) 2 is the string tension in the fun- 



damental representation. Furthermore, the temperature of the dimensional reduction T* in 
Eq. (133]) can be obtained from the equation (7f(T*) = Of, where Of(T) = [0.566^ ■ T] 2 is the 
high-temperature parametrization of the spatial string tension [13]. Solving this equation 



numerically, one gets 



15| 



1.28T,.. 



12| Plat = P\*t{T)i one 



Finally, from the lattice data for the pressure of the gluon plasma 
obtains the entropy density s = dp\ at /dT. Its plot can be found in Ref. |15| . 

With the temperature dependence of all the quantities fixed in this way, we numerically 
calculate the ratios to the entropy density of the bulk viscosities given by Eqs. (|18j) and (127]) • 
The results of this calculation are plotted in Fig. [1] For comparison, in the same Fig. [1] we 
plot the ratio Cpert/s, where the perturbative bulk viscosity 

0.443a 2 T 3 



pert 



ln(7.14/s T )" 



(34) 



with a s = g^/(4ir), was obtained in Ref. 20| in the leading logarithmic approximation. 



For illustrative purposes, we extrapolate this weak-coupling formula down to T = T c . At 



temperatures T > 2T C , our results scale as 



g% (cf. Introduction), whereas the 
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Figure 1: Calculated ratios Qt/s and °- 83 /s as functions of temperature. Also shown for com- 
parison are perturbative values Cpert/s, where Cpert is given by Eq. and extrapolated down to 
T = T C . 

perturbative result (134"|) scales as ~ g\. One observes a qualitative agreement between 
these two scaling laws and the corresponding curves in Fig. [TJ 

For the calculation of the (^ =0 ' 83 /s)-ratio, where ^£ =a83 is given by Eq. (l3~2l) . we use the 
two-loop running coupling jl3| 

^ =2&oln A + ^ 1 H Aj' Where6l = TlT6^J • 
The corresponding numerical results as a function of temperature are plotted in Fig. [2l 
The (t]t/s) -ratio, with given by Eq. (|14p . is also shown in Fig. [2] for comparison. The 
numerical values of the (t]^ =0 - S3 / s)-ratio exhibit only a small decrease due to the nonconfining 

k=0.83 

nonperturbative interactions. At temperatures T > 2T C , where ^ ~ — ~ g^, the 
calculated contributions of stochastic background fields to the shear-viscosity to the entropy- 
density ratio become subdominant compared to the contribution of valence gluons. The 
latter should gradually provide an increase of the full shear-viscosity to the entropy-density 
ratio towards the perturbative result, which is of the order of O I 4 1 ^onst I • The same 
applies to the bulk-viscosity to the entropy-density ratio, whose calculated nonperturbative 
(9(<7y)-part at T > 2T C should be gradually taken over by the perturbative O 
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Figure 2: Calculated ratio rf^ °- 83 /s as a function of temperature. Also shown for comparison are 
tjt/s, Eq. (I14p . and the conjectured lower bound of l/(47r) realized in N = 4 SYM. 

contribution of valence gluons. 



V. CONCLUDING REMARKS 



In this paper, we calculated the contributions of stochastic background fields to the shear 
and bulk viscosities of the gluon plasma in SU(3) Yang-Mills theory. These contributions 
correspond to two types of nonperturbative self-interactions of the background fields, namely 
the confining and the nonconfining ones. While the contribution of confining self-interactions 
to the shear viscosity had already been obtained in Ref. 15], here we calculated the contri- 
butions of both types of self-interactions to the both viscosities. Our method is based on 
the Kubo formulae, by means of which the correlation functions of the energy-momentum 
tensor, receiving the two above-mentioned types of nonperturbative contributions, can yield 
the spectral functions of the shear and bulk viscosities. The condition of the Matsu 
mode independence of the spectral functions, together with their Lorentzian shape 
leads to the unique correlation function of gluonic field strengths, given by Eqs. 
Remarkably, this correlation function is the same for both the bulk and the shear viscosities. 
Its amplitude, Eq. (}26l) . defines the two viscosities with the nonconfining nonperturbative 







bara- 
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self- interactions of stochastic background fields taken into account [cf. Eqs. ( |27|) and ( 132]) ]. 
Numerical results for the ratios of the viscosities to the entropy density are plotted in Figs. [1] 
and [2] They show that the found 15%-decrease of the viscosities due to the nonconfining 
nonperturbative self-interactions does not lead to somewhat significant deviations of the two 
ratios from their values corresponding to the confining self-interactions alone. The amount 
of this decrease is close to 17%, that is the relative contribution of the nonconfining part to 
the nonperturbative self-interactions of the background fields. 

At sufficiently high temperatures, T 3> T c , contributions to the viscosities produced 
by valence gluons should dominate over the above-calculated contributions of stochastic 
background fields. In particular, the purely perturbative contributions are strictly additive 
to those of the background fields 15] . Instead, at smaller temperatures, T ~ T c , a mixing 
can occur between the nonperturbative contributions of the background fields and those 
of the spatially confined valence gluons. In Ref. fl^ . such a mixing was studied for the 
pressure and the interaction measure of the gluon plasma. It was shown there that, for 
these thermodynamic quantities, spatial confinement of valence gluons plays a small role 
at temperatures T c < T < T* compared to other nonperturbative effects. For this reason, 
one can expect that its role is small for the transport coefficients as well. In particular, it 
is unlikely that the valence gluons can somewhat significantly change the obtained rapid 
decrease of the background-fields' contribution to the shear-viscosity to the entropy-density 



ratio at T c < T < T*. The change of this decrease to the O ( g4 ln CO nst ] -increase at T ^> T c 



9T 



implies the existence of a minimum of the full shear- viscosity to the entropy-density ratio at 
intermediate temperatures. Once found, the temperature at which this minimum occurs can 
be associated with a transition of the gluon plasma from the phase of a strongly interacting 
quantum liquid to the phase of a weakly interacting gas of gluons (cf. the minima of the 
shear-viscosity to the entropy-density ratio occuring nearby liquid-gas phase transitions for 
water, helium, and nitrogen, mentioned in Introduction). In the forthcoming publications, 
we plan to quantify these statements by explicitly calculating contributions of valence gluons 
to the shear and bulk viscosities. 
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